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Abstract The lithium plateau observed in halo stars has long appeared as a paradox in the general context of 
the lithium abundance behavior in stellar outer layers. First, the plateau is flat, second, the lithium abundance 
dispersion is extremely small. This seems in contradiction with the large lithium variations observed in younger 
stars. It is also difficult to understand theoretically : as lithium nuclei are destroyed by nuclear reactions at a 
relatively low temperature (= 2.5 million degrees), the occurrence of macroscopic motions in the stellar outer 
layers easily lead to lithium depletion at the surface. On the other hand, if no macroscopic motions occur in 
the stellar gas, lithium is subject to microscopic diffusion which, in the case of halo stars, should also lead to 
depletion. Several ideas have been proposed to account for the lithium behavior in halo stars. The most promising 
possi bilities were rotation al-indu ced mixin g, which could reduce lithium in the same way for all the stars (Vauclair 



1988 



Pinsonneault et al. 1992 and [19990 and mass-loss, which could oppose the lithium settling (Vauclair & 



Charbonnel 1995 and 1998). In both c ases h owever, the parameters should be tightly adjusted to prevent any 



dispersion in the final results. Vauclair (1999) (paper I) looked for a physical process which could occur in slowly 
rotating stars and explain why the dispersion of the lithium abundances in the halo stars' plateau is so small. She 
pointed out that the ^-gradient terms which appear in the computations of the meridional circulation velocity 



(e.g. Mestel 1953) were not introduced in previous computations of rotationally-induced mixing. This can lead 
to a self-regulating process which reduces the efficiency of the meridional circulation as well as the microscopic 
diffusion. Here we present numerical computations of this process and its influence on the lithium abundance 
variations in halo stars. We show that in slowly rotating stars, under some conditions, lithium can be depleted by 
a factor of up to two with a dispersion smaller than 0.1 dex in the middle part of the lithium plateau. We derive 
a primordial lithium abundance of 2.5±0.1, consistent with the recent determinations of D/H and 4 He/H. 
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1. Introduction 

Prior t o the first lithium detection in halo stars by Spite & 
Spite ( 1982 ), the primordial 7 Li value was believed to lie 
around A(Li) = 3.3 (in the log H = 12 scale), which cor- 
responds to both the meteoritic and the maximum lithium 
abundance in Pop I stars throughout the entire observed 
T e B range (8000-4500 K). Below this maximum value, the 
observations revealed a large scatter in the Li - T c g plane, 
with an evidence of lithium decrease with age for stars of 
similar T e g. 

For this reason it was expected that, due to their age, 
halo stars should have destroyed all their lithium. In this 
context, the first lithium detection in these stars, with a 
value A(Li)= 2.1, only one order of magnitude smaller 
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than the maximum Pop I value, was a surprise. The shock 
was still stronger when the observations in different stars, 
with different masses and metallicities, presented the same 
lithium value over a range of 800 K in effective tempera- 
ture (6250-5480 K). 

From then on, the observers switched to another ex- 
treme : the idea that lithium should be depleted in these 
stars was completely rejected on the basis of the constancy 
and very small dispersion of its abundances. The observed 
value was then supposed to represent the primordial abun- 
dance, without any modification, and was "sold" as such to 
cosmologists. Meanwhile, theoreticians went on claiming 
that, from all possible computations of stellar structure, 
the lithium abundance could not keep exactly the same 
value in stars during 14 billions years, without any alter- 
ation. 
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In the present paper, we argue that we can reconcile 
the theory with the observations and obtain a consistent 
view of the lithium behavior in stars. 

In section 2, we discuss recent improvements on the 
observational side and show that the careful studies of the 
lithium abundances in the " Spite plateau", obtained from 
various authors, lead to the conclusion that the dispersion 
is smaller than 0.1 dex. Meanwhile, the detection of 6 Li, if 
confirmed, puts strong constraints on the possible lithium 
destruction by nuclear reactions in these stars (otherwise 
no 6 Li should be left). 

A theoretical discussion about previous models for 
lithium abundance variations in population II stars is 
given in section 3. 

In section 4, we show how the introduction of the ef- 
fect of /^-gradients in the computati ons of rotation-induced 
mixing, as discussed by Vauclair ( 199S ) (paper I) mod- 
ifies the previously published results. We argue that a 
self-regulating process may take place, which may explain 
most of the lithium observations in halo stars. 

Section 5 is devoted to the results obtained from nu- 
merical computations of this process and its influence on 
the lithium abundance variations in different stars with 
various masses, rotational velocities and metallicities. 

We give our conclusions in section 6, including a dis- 
cussion about the assumptions introduced in our compu- 
tations. In the framework of our model, we obtain a pri- 
mordial lithium abundance of 2. 5 ±0.1. 

2. The lithium plateau : observational constraints 

Since th e first observations of the lithium plateau by Spite 
& Spite ( 1982 ), many abundance determinations have con- 
firmed the constancy and very small dispersion of the 
lithium value in most halo stars with effective tempera- 
tures larger than 5500 K. This result lead the observers to 
the conclusion that the observed lithium abundance was 
exactly the pristine value, although theoretical computa- 
tions predicted some lithium depletion in these stars. We 
first give a short summary of the important observational 
constraints and then we will discuss the theoretical models 
proposed in the literature. 

2.1. Trends with temperature and metallicity 

The observed lithium plateau is flat, at least to a 
first approximation. Some doubts still appear, however, 
about possible small variations with effective temperature 
and/or metalli city. 

Thorburn ( |1994| ) found that the lithium plateau pre- 
sented a positive slope with T c g which, she argued, would 
not exist if the lithium was primordial. She also found a 
general trend of the lithium abundance with metallicity 
which could explain part of the scatter that she observed 
on the "Spite plateau". She suggested that the slight in- 
crease of the lithium abundance with [Fe /H] could be due 
to a small production of lithium in the early life of the 
galaxy. 



While Norris et al. Ql994j ) and Ryan et al. fll996|) con- 
firmed these results, Molaro et al. ( 1995D argued that the 
slope of the plateau disappears when a subset of stars with 
temperature based on Balmer lines profiles is adopted. 
They found no correlation between lithium abundances 
and metallicity for [Fe/H] < —1.4. 

Later on, Bonifacio & Molaro ( [1997 ) studied a sam- 
ple of 41 plateau stars, to investigate in more detail the 
distribution of lithium abundances with [Fe/H] and T e g. 
They used a new effective temperature scale based on the 
new accurate infrared flux method (IRFM, Blackwell et 
al. 199C) applied to a large sample of stars by Alonso et 
al. (1996). They found only a tiny trend with T e g and no 
trend with [Fe/H]. Using mean lithium values corrected 
for standard depletion and NLTE effects, they give for 
the primordial abundance : A(Lio) — 2.238 ± 0.012. 

Meanwhile Ryan et al. ( 1999) ) studied 23 very metal- 
poor stars chosen to enable a precise measurement of the 
dispersion in the "Spite plateau" : they found no evi- 
dence for a trend with T e s but they do recover a strong 
dependence of the lithium abundance with metallicity, 
which they attribute to the chemical evolution in the early 
galaxy. The value they derive for the primordial abun- 
dance, taking into account galactic chemical evolution but 
no depletion, lies below the mean plateau value : 2.0 dex 
only, with a systematic error of 0.1 dex. 

While the question of slight lithium slopes with effec- 
tive temperature and/or metallicity is not yet completely 
settled, everyone agrees on the fact that some stars oth- 
erwise indistinguishable from normal plateau stars show 
large lithium defic iencies, well belo w the platea u itse lf 
(Hobbs et al. |l99l| , Thorburn |l994j Norris et al. |l997|) . 
Although these stars are certainly lithium depleted, so 
that their observed abundances are quite different from 
the pristine one, they must be accou nted for in the theo- 
retical scenarios. Norris et al. ( 1997 ) measured the abun- 
dances of 14 elements (Fe, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, 
Co, Ni, Sr, Y, Ba) for the metal poor stars G66-30, G139- 
8 and G186-26, which are highly lithium deficient and lie 
near the main sequence cut-off. They found no abundance 
anomalies for the other elements that one might associate 
with the lithium deficiencies. This result is in favor of a 
nuclear lithium destruction and negligible element settling 
in these stars. 



2.2. Lithium abundance dispersion in the plateau 

Since 1982, many authors have discussed the possible ex- 
istence of an intrinsic dispersion in the "Spite plateau". 



Deliyannis et al. ([1993]), supported by Thorburn ( [19940 , 
argued that the spread in measured plateau stars abun- 
dances exceeded that expected from observational data. 



They gave values o f abo ut 20% (Deliyan nis et al. 1993) 
to 25% (Thorburn |l994|) . Molaro et al ( |l995| ), Spite et 
al. (|l996|) and Bonifacio & Molaro (|l997| ) have questioned 
whether some of the error estimates in these earlier works 
were realistic and have suggested that the dispersion is 
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smaller. According to Molaro et al. ( 1995 ), all the ob- 
servations are consistent with the same pristine lithium 
abundance and the errors estimated for individual stars 



account for the observed dispersion. Spite et al. (1996) 
argued that the scatter is fully explained by the temper- 
ature equivalent width errors and that the intrinsic scat- 



ter, if real, is small. Bonifacio & Molaro (1997) revised 
the lithium abundances in a selected sample of halo stars 
using the infr ared flux method (IRFM) temperatures by 
Alonso et al. ( 1996 ). They found no evidence for intrinsic 
dispersion. 



In the sample of Ryan et al. (199S), 21 stars out of 23 
have abundances consistent with an observed spread of a 
mere 0.031 dex with reference to the metallicity trend that 
they claim to observe. Because the formal errors are 0.033 
dex, they conclude that the intrinsic spread of the Li abun- 
dances at a given metallicity is lower than 0.02 dex and 
consistent with zero at the very metal poor halo turnoff. 
This maximum dispersion is much lower than those ob- 
tained by previous studies (Spite et al. ~ 0.06-0.08 dex, 
Bonifacio and Molaro 0.07 dex) who did not introduce any 
systematic variation with metallicity. 

The situation is different in globular clusters where 
the lithium abundances in stars around the turn-off seem 
to show a dispersion l arger than that observed fo r field 
stars. Deliyannis et al. ( 1995j ) and Boesgaard et al. ( 1998 ) 
observed the stars near the turnoff of the old metal poor 
globular cluster M92. They reveal a spread in abundances 
of a factor of 2-3 for a small sample of stars. Moreover 
some of them show a high lithium abundance compared 
to the halo field plateau. 

These variations inside globular clusters may be re- 
lated to local pollution and other effects which will not be 
discussed here. In any case the very small lithium disper- 
sion in halo field stars (probably less than 10% from the 
most recent studies), acts as a strong constraint on the 
theoretical models. 

As we will see below, these results argue against strong 
lithium depletion in plateau stars. Recent observations of 
6 Li confirm this idea. 6 Li has indeed been observed in the 
atmosphere of the halo stars HP 84937 by Smith et al. 



(1993), Hobbs & Thorburn (|199j), Hobb s fc T horburn 
( |l997| ) and BD +26° 3578 by Smith et al. ( |l998| ). As this 
light lithium isotope is more easily destroyed than 7 Li, 
these observations, if confirmed, seem to prove that no 
strong lithium destruction by nuclear reactions occured in 
these stars during main sequence or pre-main sequence. 
Small destruction rates may however be allowed : we will 
see that a 7 Li depletion by a factor of two leads to a 6 Li 
decrease by less than a factor of five. 

3. The lithium plateau : theoretical discussion 

The expected lithium depletion in halo stars is primarily 
due to the diffusion processes which take place in the ra- 
diative regions inside the stars, below the outer convective 
zones. Due to pressure and temperature gradients, lithium 
settles down, as well as helium and other heavier elements. 



Mixing induced by rotation, internal waves, or mass loss 
related motions may slow down the settling process, but 
then it brings up to the convective zone matter in which 
lithium has been destroyed by nuclear reactions. 

These lithium depleting processes have been exten- 
sively studied in the literature. Michaud et al. (1984) first 
computed the lithium abundance variations in halo stars 
including element separation. Their models predicted a 
downward curvature in the Li isochrones for increasing 
effective temperature, in contradiction with the observa- 
tions. This result, simply due to the decrease of the con- 
vective depth for hotter stars, has been confi rmed many 
times in the literat ure (Deliyannis et al. 1990] , Deliyannis 
& Demarq ue |l99lj Proffitt & Michaud |l99l], C haboyer & 
Demarque 1994 , Vauclair & Charbonnel |l995|) . 

Models including rotation-induced mixing were pro- 
posed to prevent the effect of element separation. As such 
a mixing leads to lithium destruction by nuclear reactions, 
it was suggested that this destruction could be identical for 
all the sta rs, so that the lithium plat eau could be preserved 
(Vauclair |l988j, Pin sonneault et al. |l992j [1999| , Chaboyer 
& Demarque 1994 ). However, if a constant lithium abun- 
dance can be obtained in this case, it becomes extremely 
difficult to account for such a small dispersion as observed. 

Stellar models which include instabilities related to an- 
gular momentum and chemical species transport require 
additional input physics beyond standard models. The im- 
portant new ingredients include : 

1- a distribution of initial angular momenta 

2- a description for angular momentum losses 

3- a prescription for the internal transport of angular 
momentum and for the associated mixing processes in the 
stellar radiative regions. 

Us ing va rious initial angular momenta, Pinsonneault 
et al. ( 1992 ) obtained lithium depletion factors between 5 
and 10 at fixed T ff- The lithium abundances in the plateau 
slightly decreased for larger T e s in contradiction with the 
observ ations (this trend was already predicted by Vauclair 
( 1988| )). This effect however is quite attenuated for older 
ages (14 billion years) as the more lithium depleted stars 
at the hot end of the platea u mov e into the giant branch. 

Chaboyer & Demarque ( 1994 ) could also reproduce a 
lithium plateau with a large lithium depletion in a simi- 
lar framework, using different angular momentum loss and 
rotational mixing laws. They included meridional circula- 
tion, dynamical shear, secular shear, Solberg-Hoiland and 
Goldreich-Schubert-Fricke instabilities. Their results also 
showed the same trend of abundance decrease with in- 
creasing T c ff. 

Using a distribution of initial angular momenta as in- 
ferred from stellar rotation data in young open clusters, 
Pinsonneault et al. (199S) obtained a well defined, nearly 
flat lithium plateau, without the downward trend at the 
hot end. A modest scatter does remain however, which in- 
creases with the average lithium depletion : this constraint 
excludes a lithium depletion larger than a factor of 3. 

Considering the difficulty for theoretical computations, 
including element settling and rotation induced mixing, 
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to account for both the flatness of the lithium pl ateau 
and the small dispersion, Vauclair & Charbonnel (1995) 
and Vauclair & Charbonnel (199S) suggested taking into 
account a stellar wind. In some cases, a wind could pre- 
vent the element settling without bringing up to the con- 
vective zone lithium depleted matter. Computations were 
done for various mass loss rates, between 1 and 1000 times 
the solar wind. Rates larger than 10 times the solar wind 
were needed to have any effect on the lithium depletion, 
and the best results were obtained for rates of about 
10 _12 M©.?;r _1 . In this case the lithium plateau was nicely 
reproduced. For larger rates, nuclearly depleted matter be- 
gan to migrate up into the outer layers, rapidly leading to 
large lithium destruction, in contradiction with the obser- 
vations. Although appealing, this model suffers from the 
fact that it needs winds larger than the solar wind, and 
that various rates lead to a non-negligible scatter in the 
results. 

The fact that the dispersion in the lithium plateau 
is extremely small compared to the observational errors is 
the strongest constraint for the theoretical models. Several 
assumptions have lead to models which could account for 
the lithium abundance constancy, but the small observed 
scatter could only be reproduced using some "ad hoc" hy- 
pothesis on the initial parameters. 

In all these computations, however, the effect of the 
/i-gradients induced by element separation was not taken 
into account. In the present paper, we claim that when 
it is introduced in the computations of rotation-induced 
mixing, this process can explain the features observed in 
halo stars, including the very small dispersion. 



4. Diffusion and mixing in the presence of 
^-gradients 

The meridional circulation velocity in stars, in the pres- 
ence of ^-gradients, is the sum of two terms, one due to 
the classical thermal imbalance (fi-currents) and the other 
one due to the induced horizontal //-gradients (//-induced 
currents, or /i-currents in short). In the most g eneral cases, 
^-currents are opposite to 17 -currents (Mestel [l953 , Zahn 
1992| , Maeder & Zahn |1998| ). When element settling oc- 
curs below the stellar outer convective zone in cool stars, 
a small helium gradient builds, even in the presence of 
circulation. 

Then a new process must take place, which had not 
been taken into account in previous computations of 



diffusion-induced //-gradients. Mestel & Moss ( |1986[ ) gave 
a lengthy discussion of this effect for nuclearly-induced //- 
gradients. Chaboyer et al. (1995) simulated a //-gradient 
effect in their diffusion computations in the form of a re- 
ducing factor f p which was taken as a parameter. Vauclair 
( 1999| ) (paper I) showed that the resulting //- gradients are 
rapidly large enough to create //-currents of the same or- 
der as ^-currents. Then a self- regulating process may take 
place, in which both the circulation and the settling are 
strongly reduced. 



Here we give precise computations of the //-currents 
with the following simplifying assumptions : differential 
rotation is supposed negligible inside the stars (as it is 
presently in the sun) and the stellar rotation velocity is 
taken as constant along the stellar lifetime. These assump- 
tions neglect the possibility of rotational breaking in early 
stages of stellar evolution. This is discussed in the conclu- 
sion. 

4.1. Computations ofQ and fi-currents 

In the present paper, we have computed these currents in 
halo stars for different metallicities and rotation velocities 
and we have done a complete treatment of the resulting 
diffusion of the chemical species, including mixing and set- 
tling. 

As in paper I, we have neglected the deviations from 
a perfect gas law as well as the energy production terms, 
which are completely negligible in the regions of the star 
where the process takes place. 

We also assumed a negligible differential rotation, as 
observed inside the Sun from helioseismic studies. The cor- 
responding condition on Q is : 



,dlnf2, n 2 r 3 

< 



<9 In 7 



GM 



(1) 



The vertical meridional circulation velocity may be 
written (paper I) : 



ad 



Vad-V + V^ g 



(2) 



where g represents the local gravity, V a( j and V the usual 
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Here p represents the density average on the level sur- 
face (~ p) while p m is the mean density inside the sphere 
of radius r; Ht is the temperature scale height; A repre- 

M 

sents the horizontal u fluctuations — ; y M and yt represent 

M 

the derivatives : 
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While the term with the second derivative of A was ne- 
glected in paper I, these expressions are extensively com- 
puted in the present numerical computations. 

Writing u r as : 

u r = U r P 2 (cos 6) (7) 
the horizontal meridional velocity is given by : 

(8) 



ug = — — ^ ^— (pr 2 U r ) sin # cos ( 

2pr dr 



In the following, we also use the notations Un and 
defined as : 



U r = U n + Uu 



(9) 



where Un represents the part of U r which includes En, 
and U ^ the part including E^, namely : 



Un = 



Vad 



ad 



v + V, 



EnL_ 
9 M 



(10) 



idem for U^, 

The expressions and orders of magnitude of the hori- 
zontal /^-gradients A were extensively discussed in paper I. 
In all cases, on a given level surface A is proportional to 
the vertical /i-gradient. A general expression can be given 
in the form : 



A 



a \ ■ r ■ V In /x 



(11) 



where a a is a constant, which is related to the horizon- 
tal diffusion coefficient in case of shear flow instabilities by 
(Zahn 1993, paper I) : 



«A 



1 

6a h 



where 



(12) 



4.2. Element diffusion and abundance variations 



As discussed many times in the literature, stellar regions 
in which no macroscopic motions take place are subject to 
element settling induced by pressure and te mper ature gra- 
dients and radiative acceleration ( Michaud 197C , Michaud 
et al. |1976| (M CV2) Vauc lair & Vauclair |l982j Vauclair 



& Charbonnel 1995| , 1998 ). When macroscopic motions 
occur, they slow down the settling but do not prevent 
it completely until an equi librium concen tration gradient 
is reached (V auclai r et al. |l978a| (V2M), |l978b| (V2SM), 
Richard et al. |1996| ). 

In paper I we showed that in halo stars this equilibrium 
gradient is never reached before the ^-currents become 
of the same order of magnitude as the f2-currents. Thus 
the abundance variations induced by settling and rotation- 
induced mixing must always be computed simultaneously. 

Computations of element settling have been described 
many times. He re we use the same prescriptions as in 
Richard et al. (11996D with collision cross sections and 



al. ( [1986 ). Radiative accelerations are neglected, as usual 
in these cool stars (precise computations in the Sun by 
Turcotte et al. ( 1998| ) show that they are indeed negligi- 
ble for our purpose). 

Rotation-induced mixing is also computed as in 
Richard et al. ( 1996| ) except for the introduction of 
the /^-currents in the velocity of meridional circulation. 
Following Zahn (1992) and ( 1993[ ), meridional circulation 
is suppose to induce shear flow instabilities which, in a 
density stratified medium, lead to large horizontal diffu- 
sivities. The combination of the circulation and horizontal 
mixing may be treated as a vertical effective diffusion pro- 
cess, described with an effective diffusion coefficient : 



D cS 



C h 



U r 



with C h ~ 30 



(13) 



Abundance variations of helium, lithium and heavier 
elements are computed within this framework, all along 
evolutionary sequences. The modifications of the internal 
structure induced by element diffusion are taken into ac- 



count in the computations, as in Richard et al. (1996) 



microscopic diffusion coefficients as given by Paquette et 



4.3. Self-regulating process 

The originality of the present paper lies in the introduction 
of /^-currents in the meridional circulation. As discussed 
in paper I and verified in the complete numerical compu- 
tations (see below), /i-currents become of the same order 
as ^-currents for small ^-gradients, rapidly reached from 
helium settling below the convective zone. In this case we 
claim that, under some conditions, a self-regulating pro- 
cess may take place in which both the circulation and the 
settling are strongly reduced. 

At the beginning of the stellar evolution, |E M is smaller 
than | -En | ; meridional circulation proceeds and induces 
transport of chemical elements coupled with settling. The 
helium concentration gradient below the convective zone 
increases (in absolute value) leading to an increasing ver- 
tical ^-gradient . The horizontal /i-gradient A follows pro- 
portionally until \Ep\ becomes of the same order as |£n|, 
where the circulation velocity is expected to vanish. 

Then the whole process is modified and does not follow 
any longer the classical treatment of the Eddington-Sweet 
circulation. Let us call A cr u the horizontal /^-gradient for 
which \E^\ = | -En |- Without microscopic diffusion, the 
whole circulation would become frozen, keeping horizontal 
/^-gradients of the order of A cr it in every layer. However 
microscopic diffusion still proceeds below the convective 
zone, leading to a slight decrease of A in the upper radia- 
tive layers (this is due to the fact that there is no horizon- 
tal /x-gradient inside the convective zone because of rapid 
mixing) . 

As A becomes smaller than A crit , circulation proceeds 
just below the convective zone, in the direction of Uq. 
Doing so, due to mass conservation, it lifts up matter from 
below in an ascending flow and pushes it down in a de- 
scending flow. Although this description looks like that 
of normal meridional circulation, there is a fundamental 
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difference : in the deeper layers, where A was close to 
A cr jt , this extra motion creates an over-critical horizontal 
//-gradient, which leads to a local horizontal motion in the 
direction of U^. Then part of the matter which falls down 
in the descending flow is transferred to the ascending flow 
and goes back up to the convective zone where it came 
from. We can infer that, when this occurs, the global el- 
ement depletion from the convective zone is reduced in 
a significant way, leading to time scales larger than the 
stellar ages. 

A good treatment of this self-regulating process needs 
a complete numerical simulation, which is planned for a 
forthcoming paper. Here we used the following approxi- 
mation. 

In the layers where the circulation time scale (that is, 
the time scale of both the ^-currents and the //-currents 
taken separately) is smaller than the settling time scale, 
we assumed that, as soon as [E^[ becomes of the same 
order as \Eq\ (within 10%), the //-gradient remains equal 
to the equilibrium one. 

On the contrary, in the layers where the circulation 
time scale is larger than the settling time scale, we as- 
sumed that the readjustement would have no time to take 
place in a settling time scale and element depletion pro- 
ceeds unaltered. 

In other words, the //-gradient remains constant below 
the convective zone, where Un is larger than the settling 
velocity, as soon as \E^\ ~ \Eq\. Meanwhile the overall 
helium and lithium abundances are still subject to deple- 
tion due to diffusion deeper in the star, at the place where 
the settling velocity becomes larger than Uq (see below, 
Figure 2). 



5. Computational results 

5.1. Evolution of fi-currents with time in low 
metallicity stars 

First we present the results obtained for four low metallic- 
ity stars with masses of 0.75Mq (hot end of the plateau), 
0.70M© (middle of the plateau) ,0.65M© (cool end of the 
plateau) and 0.60M Q to show the behavior of the //- 
currents with time. Here the rotation velocity is assumed 
constant, equal to V ro t = 5 km.s -1 and the metallicity is 
chosen equal to [Fe/H]=-2.0. 

Figures 1 (a, b, c, d) display the \Eq\ and [E^[ pro- 
files below the convective zones at different evolutionary 
stages. Each star arrives on the main sequence with nearly 
homogenous composition. At that time, \E^\ is smaller 
than \Eq\ everywhere inside the star. Meridional circula- 
tion and element settling can take place in the radiative 
regions of the star as described above. As a consequence 
of the helium settling, the //-currents increase below the 
convective zone and become rapidly of the same order as 
ri-currents. 



5.2. Lithium abundance variations in plateau stars 

The lithium abundance variations with time have been 
computed using the following scenario : microscopic dif- 
fusion damped by rotation-induced mixing occurs nor- 
mally until the relative difference between \E^\ and \Eq\ 
becomes smaller than = 10%. Then the self- regulating 
process is supposed to take place, in which diffusion and 
mixing are strongly reduced (with time scales larger than 
the stellar lifetime). This occurs down to the place where 
the microscopic diffusion velocity becomes preponderant 
(Figure 2, a, b, c). Below this layer, diffusion is assumed 
to proceed freely. 

Figure 3 presents the lithium abundance variations 
with time for a bunch of different rotational velocities in 
a 0.70M Q star. For small rotation velocities, the lithium 
abundance variations are mainly due to microscopic dif- 
fusion, decreasing for increasing rotation velocity. For ve- 
locities larger than 6 km.s -1 the depletion is mainly due 
to nuclear destruction, which now increases with increas- 
ing velocity. For velocities between 2.5 and 7 km.s -1 , the 
resulting dispersion is of order 0.06 dex, going up to .12 
dex if the 7.5 km.s -1 possibility is included. 

Figure 4 (a, b, c) presents the lithium abundance 
variations with time in the outer layers of the four low- 
metallicity stars ([Fe/H] = -2) with three different rota- 
tion velocities: 2.5, 5.0 and 7.5 km.s -1 . The 0.65M Q and 
0.70M Q stars show very similar abundances and their vari- 
ations with the rotation velocity is small except for the 
largest value of 7.5 km.s -1 . In the 0.75M Q star the lithium 
abundance is more sensitive to the rotation velocity due 
to the smaller diffusion time scale. For small velocities the 
lithium value in the O.75M star is smaller than in cooler 
stars while it is larger for velocities larger than 5 km.s -1 . 
For the cooler 0.60M Q star, lithium is always more de- 
pleted due to the small distance between the bottom of 
the convective zone and the nuclear destruction layer. 

Figure 5 (a, b, c, d) displays, as an example, the abun- 
dance profiles for lithium inside the four halo stars, for a 
rotation velocity of 6 km.s -1 . 

5.3. Influence of metallicity and mixing length 
parameter on the lithium abundances 

We have tested the influence of metallicity on the lithium 
abundance variations for all the stars, with two different 
metallicities [Fe/H] = -2.0 and [Fe/H] = -3.5. The 
curves representing the lithium variations with time are 
plotted in Figure 6. These curves are very close, which 
leads to the conclusion of a very small dependence of the 
lithium abundance on metallicity. 

We have also tested the influence of varying the mixing 
length parameter used for the computation of the convec- 
tive zone: the results are presented in Figure 7. The differ- 
ences are very small : the results are not sensitive to small 
variations of the depth of the convective zone. 
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5.4. The case of cool stars 
tachocline 



introduction of a 



The lithium depletion obtained for the O.6OM star, which 
is cooler than the lithium plateau, is not sufficient to ac- 
count for the observations, except for a rotation velocity 
of 7.5 km.s -1 . The observations show, for these effective 
temperatures, a lithium depletion by at least a factor of 
10, while the computations obtained by using the same 
physical inputs as those introduced for the three others 
stars give a lithium depletion of only 60% for 2.5 km.s" 1 
and 80% for 5 km.s" 1 . 

However, at such effective temperatures, a tachocline 
layer is supposed to take place below the convective zone, 



as in the Sun (Spiegel & Zahn 1992). We introduced in the 
computations an exponential mixing just below the con- 
vective region to reproduce the effects of this tachocline. 
The width of the mixed layer was taken as 0.01 R* and the 
diffusion coefficient was 100 cm 2 .s -1 below the convective 
zone, decreasing exponentially underneath. 

We then obtained the lithium variation with time as 
plotted on Figure 8 for a velocity of 5 km.s" 1 , giving a 
lithium depletion by a factor 20, consistent with the ob- 
servations. We have further checked that such a narrow 
tachocline layer, if present in hotter stars, would not mod- 
ify the previously obtained results. 



5.5. Abundance variations of 6 Li 

We have done the same computations for 6 Li as for 7 Li 
in the plateau stars. Figure 9 displays the results for a 
0.70M© star with metallicity [Fe/H] = -2 and the same 
values of the rotation velocity as for 7 Li in Figure 3. The 
results differ for the two isotopes because 6 Li is more easily 
destroyed by nuclear reactions than 7 Li, while the influ- 
ence of microscopic diffusion is about the same. For small 
rotation velocities, less than 7 km.s -1 , 6 Li is destroyed by 
about a factor of two. For larger velocities it goes up to a 
factor of 5. In any case, if 6 Li is present at the beginning 
of the stellar evolution, a non negligible amount remains 
at the age of halo stars. 



6. Discussion 

The results obtained in this paper are compared to the 
observat ional points and error bars given by Bonifacio & 
Molaro ( [1997D in Figure 10 (a) and (b). The solid curve 
drawn in Figure 10(a) represents averaged depletion values 
for rotation velocities between 2.5 and 7.5 km.s -1 . The 
original lithium value is A(Li) = 2.5. In figure 10 (b), 
two extreme curves are given, again with the assumption 
of rotation velocities between 2.5 and 7.5 km.s -1 . The 
expected dispersion is very small (of order 0.1 dex) at 
the middle of the lithium plateau. It becomes larger at 
both extremes: for the hotter side, it is due to the shorter 
diffusion time scale and for the cooler side to the proximity 
of the nuclear destruction layers. 



A complete treatment of the self-regulating process de- 
scribed in section 4-3 would need numerical simulations 
which are in preparation. Here we have assumed that this 
effect significantly reduces both the circulation and the 
diffusion as the circulation time scale is smaller than the 
diffusion time scale while it is supposed to proceed unal- 
tered below. Such a simple scheme gives results consistent 
with the lithium observations in halo stars. No mass loss 
has been taken into account, which means that the mass 
loss rate is assumed smaller than lO -13 M .yr -1 (Vauclair 



& Charbonnel |l995[) . 

Two important assumptions have been used in these 
computations : the stellar rotation velocity was taken as 
a constant during evolution and the differential rotation 
inside the stars was assumed negligible. 

The first assumption is impossible to check directly, as 
we have no access to the past rotational history of halo 
stars. It is generally assumed that they suffer the same 
rotational breaking as pop I stars. However they have not 
been formed in the same way, not being in the same galac- 
tic sites, and their kinematics is different. There is no rea- 
son why halo stars should behave exactly in the same way 
as open cluster stars. If the rotational velocity of halo stars 
was larger in the past, we would expect a larger lithium 
destruction before the self-regulating process we have de- 
scribed takes place. This cannot be excluded, so that the 
result we give below has to be taken as the lower limit of 
the primordial value. 

The second assumption is supported by recent helio- 
seismic results, which show that the internal rotation of 
the Sun follows a solid body law except in the convec- 
tive zone. We may suppose that the same occurs in halo 
stars. The fact that the obtained results nicely reproduce 
the observations argue in favour of this hypothesis. Note 
that the most efficient way to reach this situation seems 
to be the p resence of a very small magne tic fie ld (Mestel 
et al. 1987 , Charbo nneau & MacGregor 1992 ) which, as 
discussed by Mestel ( 1999 ), would suppress differential ro- 
tation without altering large-scale motion. Such a physics 
has to be further studied. 

It is possible, however, that in some specific stars, due 
to gravitational interaction with other stars or to some un- 
known other process, differential rotation does occur. This 
could explain the existence of some stars with no observed 
lithium below the plateau : lithium would have been de- 
stroyed there by nuclear reactions due to extra mixing. 
It could also explain why some dispersion is observed in 
globular cluster stars whil e it is not observed in field stars 
(Boesgaard & Deliyannis 1998} ) . 

Futher computations have to be done to describe more 
precisely the hydrodynamical situation which occurs when 
the ^-currents oppose the f2-currents below the convective 
zones of slowly rotating stars. Applications to lithium in 
open clusters are underway. The results presented here 
are very encouraging in this respect. They show that the 
plateau can be reproduced with no "ad hoc" adjustment of 
any parameter, with an overall lithium depletion slightly 
smaller than a factor two. In this case 6 Li is still present 
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Masse 


M=0.75 M 


M=0.70 M 


M=0.65 M 


M=0.60 M 


Li 

log yt- 

Lm 




Li 

logy— 

Lin 




Li 

logy^ 

Lm 




Li 

log YT 

Lin 




V=2.5 km.s" 1 


10 Gyrs 


-0.29 


6393 


-0.20 


6059 


-0.22 


5718 


-0.52 


5383 


12 Gyrs 


-0.38 


6524 


-0.26 


6139 


-0.28 


5773 


-0.57 


5421 


V=5.0 km.s" 1 


10 Gyrs 


-0.26 


6391 


-0.18 


6059 


-0.18 


5715 


-0.64 


5185 


12 Gyrs 


-0.31 


6473 


-0.22 


6133 


-0.21 


5770 


-0.73 


5226 


V=7.5 km.s" 1 


10 Gyrs 


-0.21 


6397 


-0.28 


6059 


-0.46 


5720 


-1.34 


5385 


12 Gyrs 


-0.23 


6507 


-0.32 


6136 


-0.52 


5774 


-1.47 


5423 



Table 1. Lithium depletion and effective temperatures of stars for the four considered masses, three different rotation 
velocities and two ages. Here no tachocline is included in the models. 



in the plateau stars. The primordial lithium value is then 
found as 2.5 ± 0.1, consistent with the recent determi- 
nations of D/H and 4 He/H (see Buries & Tytler |l99g|) . 
These results give a baryonic number r] = (5 ± 1).10~ 10 
and a baryonic density VLbh 2 between 0.015 and 0.022. 
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Figure 1. a. Evolution of ^-currents with time inside the star. The graphs show the variations with depth of both 
| £h | and {E^ (defined in eq. 4 and 5) in a 0.75M Q halo star with [Fe/H] = —2 and V rot — 5 km.s -1 . 
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Figure 2. a. Microscopic diffusion velocity below the convective zone of a 0.75M Q star compared to the classical 
meridional circulation velocity Un (defined in eq. 10). When Vdif < Un, the self- regulating process can occur while 
when Vdif > Un, the circulation cannot compensate for the effect of diffusion. 
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Age 

Figure 3. Lithium abundance variations at the surface of a 0.70M© star for rotation velocities between 2.5 and 7.5 
km.s" 1 . The minimum depletion is obtained for 6 km.s . For smaller velocities, depletion is mainly due to microscopic 
diffusion while for larger velocities it is mainly due to nuclear destruction. Due to the self-regulating process, the 
obtained dispersion is quite small. 
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Figure 4. Lithium abundance variations at the surface of the four halo stars, with three different rotation velocities 
(a) : 2.5 km.s" 1 , (b) : 5 km.s" 1 , (c) : 7.5 km-s" 1 . 
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Figure 5. Lithium profiles inside the four halo stars, with a rotation velocity of 6 km.s . Dashed line : 0.029 Gyrs, 
long dashed dotted line : 1.869 Gyrs, dotted line : 3.738 Gyrs and solid line : 10.095 Gyrs. 
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Figure 7. Lithium abundance variation with time for 0.70M Q halo stars with two different values of the mixing length 
parameter. 
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Figure 8. Lithium abundance variation with time for 0.60M Q halo stars with and without tachocline. Solid line : 
stellar model without tachocline, dashed line : stellar model with tachocline. 




Figure 9. Lithium 6 abundance variations at the surface of a 0.70M Q star for rotation velocities between 2.5 and 7.5 
km.s" 1 . 
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Figure 10. C omparison between the computations of lithium depletion and the observational values of Bonifacio & 
Molaro ( 1997 ). (a) : average curve obtained for rotation velocities between 2.5 and 7.5 km.s -1 . Here a tachocline is 
introduced as explained in section 5.4. The original lithium value is choosen as 2.5 to fit the observations, (b) : expected 
dispersion if the rotation velocities lie between 2.5 and 7.5 km.s -1 . The observational error bars are not shown for 
clarity. 



